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NATIONAL ADVISORY- COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2213 

AERODYNAMIC COEFFICIENTS FOR AN OSCILLATING AIRFOIL WITH 
HINGED FLAP, WITH TABLES FOR A MACH NUMBER OF 0.7 
By M. J. Turner and S. Rabincrwitz 

SUMMARY 


Dietze's method for the solution of Possio' s integral equation has 
been used to determine the<> chordwise distribution of lift on an oscil- 
lating airfoil with simple hinged flap in two-dimensional compressible 
flow ( subsonic) . The results of these calculations have been used to 
prepare tables of aerodynamic coefficients for lift, pitching moment 
(referred to quarter-chord point), and flap hinge moment for a Mach 
number of 0 . 7 ; the motions considered are vertical translation, airfoil 
rotation about the quarter-chord point, and rotation of the flap about 
its hinge line. 

Aerodynamic coefficients are tabulated for 4 values of t b (ratio 
of flap chord to total chord) and for 12 values of reduced frequency ojy, 
covering the range from 0 to 0.7. Results are given for one value of 
Mach number, M = 0.7. Data of this kind have already been presented 
by Dietze for one value of -the ratio of flap chord to total chord, 

Tg = 0 . 15 ; these results have been checked independently, and calcula- 
tions have been carried out for three additional values, t b = 0.24, 

O. 33 , and 0.42. Certain auxiliary parameters, which will be needed in 
any further calculations of this type, are presented for future reference. 


INTRODUCTION 


The fundamental integral equation for the pressure distribution on 
an oscillating thin airfoil moving at subsonic speed has been derived 
by Possio in reference 1. Collocation procedures have been used by 
Possio, Frazer and Skan, and others to obtain lift and moment on an 
oscillating flat plate. An important contribution has been made by 
Dietze (see references 2 and 3), who has developed an iterative procedure 
for numerical solution of Possio' s integral equation. This procedure is 
particularly well adapted to the calculation of aerodynamic loading on 
an osc ill at ing airfoil with hinged flap. 
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It has "been pointed out correctly "by Karp and Weil ( reference 
p. 11) that Dietze's procedure does not properly account for the loga- 
rithmic singularity in pressure distribution at the flap hinge-. However, 
for applications in flutter analysis the principal objective is to deter- 
mine resultant lift and moments rather than pressure distribution. Mathe- 
matically exact solutions in closed foim are available for the stationary 
thin airfoil with deflected flap, and it is found that lift and moments 
obtained by Dietze * s iterative procedure are in excellent agreement with 
theoretically exact values . . There is good reason to expect that equally 
satisfactory results will be obtained for the airfoil with oscillating 
flap. The existence of the singularity in pressure distribution is a 
consequence of the sharp comer in the idealized, broken-line profile, 
and of the associated discontinuity in downwash velocity at the flap 
hinge. The introduction of a f ini te cosine series for the downwash is 
in effect equivalent to a slight modification of the idealized profile 
by roun ding off the comer at the flap hinge. ° 

This work has been performed at Chance V ought Aircraft, under the 
sponsorship of the Bureau of Aeronautics, Navy Department, in order to 
provide data for the calculation of compressibility effects in control- 
surface flutter problems. It has been made available to the National 
Advisory C ommi ttee for Aeronautics for publication because of its general 
interest. 


SYMBOLS 


T R 

ratio of flap chord length to total chord length 

*a 

airfoil region in x,z-plane 

8y(x,t) 

vertical displacement of point on idealized profile, 
positive upward 

8y(!) 

nondimens ional representation of instantaneous chordwise 
distribution of vertical displacement 

1 

total chord length . 

0) 

0 circular frequency , 

t 

time 


x 


chordwise coordinate 
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y 

z 

5 

Vy(x,t) 

g(|) 

V 

7(5) 

7±nc 

K(s,M) 

u,v 

s 

M 

\i = 1 - \Ji 

p 

a(x,t) 

*(®r) 

AK(s,M) 

AK x (s,M) 

AK2(s,M) 

k ij 

k2n 


vertical coordinate 
spanwise coordinate 

dimensionless chordwise coordinate (2 x/l) 

vertical component of fluid velocity adjacent to airfoil 

function representing instantaneous distribution of 
vertical fluid velocity adjacent to airfoil 

(vy(x,t) = gC^e 1 ^) 

reduced frequency (o>z/2V) 
velocity of flight 

function representing distribution of dipole l ine s 

/ 

distribution of dipole lines for incompressible flow 

kernel of Possio’s integral equation 
variables of integration 
auxiliary variable ^ao r ( £ - 5o)) 

Mach number 

- M 2 

air density 
lift per unit area 

function defined by Kussner and Schwarz 

kernel difference (k(s,M) - K(s,0)) 

singular part of kernel difference 

nonsingular part of kernel difference 

constant occurring in formula for AK^ 

coefficient in polynominal representation for AK2 
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Pile* € nv 

AP & 

AMd 

AMr 

% 

% 

9R 

c gh> k gh 


coefficients in recursion formulas for solution of 
Possio 's integral equation 

lift force on airfoil strip of width. Az 

pitching moment on airfoil strip of width Az 

hinge mome nt on flap for strip of width Az 

downward displacement of quarter-chord point divided by 
semichord 

rotation of airfoil, positive in stalling direction 
rotation of flap, positive in stalling direction 
aerodynamic coefficients, where g, h=S, D, R 


BASIC THEORY 


A very complete digest of the literature on aerodynamic theory of 
oscillating airfoils has been presented by Karp, Shu, and Weil in refer- 
ence 5* Consequently, the basic theory is outlined briefly merely to 
exhibit the essential features of the computational scheme and to point 
out certain errors which have been discovered in Dietze*s formulas. 

The usual assumptions of thin-airfoil theory are adopted, leading 
to Possio* s integral equation for the chordwise lift distribution on 
the oscillating airfoil. Rectangular coordinates are employed, with 
the x-ax±s alined in the direction of the undisturbed flow. The airfoil 
is replaced by a deformable sheet of zero thickness which, in its undis- 
turbed position, occupies the region Hq^ 

- 1 / 2 = x = Z/2, y = 0 


of the x,z-plane. 

lifting surface’ executes sinusoidal oscillations in which each 
point moves along a line parallel to the vertical y-axis. Displacements 
are independent of the spanwise coordinate z, and may be represented in 
the form 


icut 


B y = B y(*,t) = ^ 5y(i)e' 
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In accordance with the usual convention, it is the real part of equa- 
tion (l) which has physical significance.. 

The y-component of velocity of a fluid particle adjacent to the 
lifting surface is related to 8 y "by the equation 


^By 2V ^5y 

v y = ^ 2 + v^ = 5 ^ + T3! 2 


( 2 ) 


where 5 = 2x/l and - 1 = | = 1 inside Rg. From equations (l) and 
( 2 ) it follows that 

y y = gtoe*"* 


(3) 


where 


sU) '= v(l%S y + ^) 


(4) 


where 


ojj. = ai2/2 V 


In Dietze's derivation of the basic integral equation the region R a 
is covered with dipole lines of density V 7 (x)e iu:ft per unit length in 
the chordwise direction. By equating the vertical velocity induced by 
the dipole covering to that given by equation ( 3 ) the following integral 
equation is obtained for the determination of 7 : 

g(|) =<D r 'j) 7 (| 0 )k(s,M) d £ 0 (5) 

subject to the condition that 7 ( 1 ) shall be finite; the kernel of the 
integral equation is given by 


K(s,M) =- 


,isXM 


k\Jl - M 2 


|h 0 ( 2 )(|s|x) - m ~~ %( 2 )( |S|X) - 

1 8 I 


. X 

i(l - #)e Ti 


,S H 


log. 


M 


Tt \jl - M 2 1 - V X - M 2 


eiu Ho (2)( (u ,j I ) 

du 


( 6 ) 



6 


MCA TN 2213 


with 


s = 0^(1 - |o) 

X = M/(l - MS) 

The lift a per unit area (positive upward) is given hy 

a(x,t) = pVrCxJe 3 ^ 


(7) 


NUMERICAL SOLUTION OF POSSIO’S INTEGRAL EQUATION 


Dietze’s approximate solution of equation (5) (Possio’s integral 
equation) is of the form 


7 * 7^ + 7i + 7 2 + 


+ 7 r 


( 8 ) 


where 7^ and 7 V (v = 1, 2, . 
equations 

r 1 


. n) are solutions of the integral 


g(5) = a* y 7 inc(So)K(s,0) d£ 0 

P 1 


(9) 


Sy(S) — ^ ) 7y(^o)K(s,0) dig* V — 1; 2; . . • * n (lO) 

J-l 


and where 


nl 


gl(l) -a**) 7inc(5o)|K( s ^ 0 ) - K(s,M)] <*t 0 

J-l 


( 11 ) 


il 

gyU) = “r^> 7 v _i(lo)[ K ( B #°) - K(s,M)] d|o, v = 2, 3, • • n 

J-l (12) 
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K(s,0) = 11m K(s,M) 
M — > 0 


The required solutions of equations (9) and (10) are obtained by the 
methods of reference 6. Convergence of the process has been proved only 
for the stationary case. However, computational experience furnishes 
convincing evidence that the process does converge in the more general 
case. 


It will be observed that Dietze's process requires essentially the . 
solution of a succession of integral equations with kernel K(s,0) from 
the incompressible problem. The functions g v (|) are obtained by 

direct integration in accordance with equations (ll) and (12) . 

In case gy(|) can be represented by a cosine series of form 

gyd) = v(Aq + 2 X An cos n 0 ^ ^ it ( 13 ) 

with 


| = -cos 0 


then it is known from the work of Kussner and Schwarz (reference 6) that 


where 


7 V U) 




®0 ■ ( ~ ~ 2 ~ ~ ) ( A 0 “ Al ) + A 1 

®n = r~~( A n-l " An+l) “ A n> n = 1 


(14) 


(15) 


and Tfcoj.) is the function of reduced frequency defined in reference 6. 
The T-function is related to Theodorsen*s C-f unction by the equa- 
tion T = 2G - 1. 


In order to facilitate the evaluation of the integral occurring in 
equation (12) the kernel difference 


AK(s,M) = K(s,M) - K(s,0) 
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is expressed in the form 

AK(s,M) = ZtfC^(s,M) + £K^(s,M) ’ ( x 6) 

where 


AKi(s,M) = _ + kll + kl2 loge j g{ , + s ^ i3 + klk lQge j B |J (1?) 
k l° - ~(l - \[l- M2) 

krL ' "*(vi - M2 " X ) " nap ‘ l0& 2(l ? - M2) ] + 


k 12 = 


^3~S 


|-1 + loge 


7M 


1 + \/l - MS (l _ M2)3/ 2 

§m* - (a . § M a) loge .. i 

J L 


H 


M2 - 


3M 2 - 2 
2(1 - M2) 3/2 



3M 2 - 2 

2(1 - MS) 3/2 


l°Se^ — 0*57722 (Euler’s constant) 
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The nonsingular part AKgCs^M) is replaced hy a polynomial of ninth 
degree-*- 

AK2(s,M) js - kpnS n ( 18 ) 

n=2 


whose coefficients are determined separately for each Mach number by 
fitting the polynomial to tabulated values in the interval |s| = 1.8. 


Equation (12) may be written in the - form 


gyU) = 7 V _! AK(s,M) d£ 0 

(19) 

If it be assumed that 


7 V _ X = “2V^p 0 cot | + 2 p n sin n^ 

(20) 

then it follows that (upon carrying out the required integrals 
defining g v in accordance with equations ( 17 ), ( 18 ), and ( 19 ) 
applying equations ( 15 ) to solve equation (10)), 

and 

7 V = ~ 2V (%o cot | + 2 £ q n sin n^ 

(21) 

where 

^n = ^ln + 92n 

(22) 


^Dietze's definition of the approximating polynomial (table 3 o f 
reference 2) differs in sign from that given here. However, it has been 
found by carefully checking the derivations that the minus sign is 
required in equation (l8) in order to justify the recursion formulas 
for calculating 7 V from 7 v _i. Tt is believed that this is merely an 

error of presentation, since Dietze's numerical results are found to be 
substantially correct. 1 



10 


MCA TN 2213 


and 


I 

»H> “ ** (-®1 + PooV + f>01 p l' + P 02 P 2 ') 

«u - K p i + “uPo' + P ii p i’ + h&z' + f>i3V’) 

q 12 = p ( p 2 + ^O-Po' + ^21 p 2* + ^22^2 ' ' + %3 P 3") 

'bn = P( P n + Vn' + V." + Y«' ")’ “ ^ 3 

ho = 2lt p v' s 0V 

10-n 

ha ~ ^ _1 ^ 1 ®*^jj- p v ,, nv' n = 

with 

p o' - -\ ( p 0 + p l) 
p l' ”-|( p 0 +p 2) 

4 

V - (Pn-l - p n+l) J n " 2 
p n" *s( p n-l' - p n+l')’ “ ® 2 

V "( p n-l" - p n + l"> p -3 


(23) 


TTip quantities a^, p^, and are defined in the translation 

of reference 2 (see table 5> P* 28, and tables 6 and 7* PP» 30-32). It 
was found that the formula for is given incorrectly by Dietze; it 

e 0 V = * T ) (Sqv + 5lv) - SlV 


should read 


m 
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Presumably this is merely an error in presentation, since, as already 
noted, Dietze * s numerical results have been verified by independent 
calculations. Numerical values of the parameters o^, ft-;> for M = 0.7 
and reduced frequencies ranging from 0 to 0.7 are given in table I. 


EVALUATION OF KERNEL AND KERNEL DIFFERENCE 


By making use of the relations among Bessel, Haakel, and Neumann 
functions and by separating the kernel into real and imaginary parts an 
expression of the following form is obtained: 


K(s,M) = K‘(s,M) + ££*'(8,14) 


-A(s,M) cos (bAM) + B(b,M) Bin (sXM) 
k\Jl - M 2 


i B(s,M) cos (sXM) - A(s,M) sin (sXM) 

k\Jl - M 2 


(25) 


where 


A = J 0 ( | s |1) - mJU N-jUsIX) - |\/l - M 2 l 0ge sin (s - 

1 * 1 - VI - M 2 7 


(l - M 2 ) sin ^s »J„ [cos u Jo(|u|M) + sin u Nq(|u|M)J du - 

(l - M?) cos -^s ^ ^ [cos u Nq(|u|M) - sin u Jq( |u|M.)]J du 


B =* N 0 (| s|X) + M JlOsIX) + |\/l - M 2 loge cos (s 

(l - M 2 ) cos (■ » r [cos u Jq(|u|M) + sin u No(|u|M)J du - 


(26) 


> s™ 
M 


(l - m 2) sin ^s ^ | Qcos u No(|u|M) - sin u Jo( |u|M)| du 
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By Introduc ing a new variable of integration 


M 

v = x u 


equations ( 2 6 ) are transformed into 
A = Jo(lslX) - M^-Ni(|s|>.) - |\/l - M 2 logg M 


ir 


1 - \ll - M 2 


=sr sin 


( s h)- 


sin 


(■ jjj) J [cos (v J 0 (|v|X) + sin (v M*)J 


dv - 


JO 

ns 


cos 


(s ^ [os (v ^B 0 (|v|X) - sin ^v |^J 0 ( |v|X) dv ( 27 ) 


B = N 0 (|s|X) + M J]_( |s|A.) + §\/l - M 2 loge 


M 


cos 


1 - \/l - M 2 

IS 

( s s) [ os ( T If) J °( l T l x ) + sin ( T ,V,X )J dT - 


(■>) 


sin 


( s i) I cos ( T I) N °^ i v i x ^ 1 sin ( v m) j °^ ,v|X) dv 


( 28 ) 



MCA TN 2213 . 


13 


Numerical values of the following integrals are required: 

I x (s) 

I 2 (s) 

I 3 (s) 

Ilt(o) 

The evaluation of 1^ and Ig can he obtained by numerical integration 
in a straightforward manner. However, since No(x) has a logarithmic 
singularity at x = 0, it is necessary to express 1 3 and 1]^ in 
different form. Upon making the substitution (reference 7, pp/ 130, 132) 

\ N 0 (x) = J 0 (x) logg H _ Bq(x) 




where 
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and integrating "by parts, 
| I 3 (s) = I x (s) logg ~ 

\ I^Cs) = 3^(s) logg 


the following equations 


Il(v) 


ns 


dv - 


cps 


UO 


ns 


^2^ t ) 


OS 


dv - 


sin 


JO » 


are obtained: 

( v m) b o( vX ) dv 

(y^)Bo(v\) dv 


(30) 


(31) 


It follows from equations (29) that Ii and I3 aria odd functions, 
while I 2 and I4 are even. The integrals occurring in equations (30) 

and ( 3l) can he evaluated numerically without difficulty; it should he 
noted that 


1 


0 


In computing numerical values of the kernel Dietze has used the 
tables of Bessel and Neumann functions given in reference 7, which do 
not permit a satisfactory determ.1 nation of Bo since No is tabulated 
to only four (in some cases three) places. In recalculating the kernel 
the seven-plane tables given in reference 8 have been used. 

For evaluation of the .integrals I n the formulas given in refer- 
ence 9, page 227, have been extended to include fifth differences. These 
formulas have been used with an interval A (sX) = 0.05, and the results 
have been checked Tip to sX = 0.30 by using an interval A(sX) = 0.02. 
Becalculated values of K(s,0.7), K(s,0), AK(s,0.7), AKi(s,0.7), and 
AKg(s,0.7) are given in table H. These values are found to be in 
close agreement with those given by Dietze; where differences exist, the 
new values are believed to be more accurate. In making comparisons it 
should be noted that Dietze has tabulated -K(s,M). 

The coefficients are obtained by fitting a ninth-degree 

polynomial (see equation (18) ) to the tabulated values of AKp in the 


11m 
v— >0 

lim 


Il(v) 


I2( T ) 
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interval |s| = 1.8 by the method of least squares.. The values obtained 
in this way for M = 0.7 axe 

*22 = -0.01*6728 - 0. 0^597^1 
N 3 = 0.023019 - 0. 023^911 
*24 = 0.009818 + 0.052855i 
k 2 5 - -0.020228 + 0. 003l|88i 

kgg = -0.00101*0 - 0.0215101 

kgY = 0.007272 - 0.0002831 
kgg = 0.000053 + 0.0031171 
N9 = -0.000997 + 0.0000121 


Since |s| = 0^(5 - £q) U " lo | = 2 the approximate representa- 

tion for AK 2 is valid' for reduced frequencies in the range 0 = aoj. = 0.9. 
Although the values of the coefficients, k p^ given here differ from those 
given by Dietze (see translation of reference 2, p. 26), there is reason- 
able agreement of coefficients for lower powers of s. Also the algebraic 
signs agree; this gives further indication that Dietze must have intro- 
duced a minus sign in his definition of the approximate representation 
for AN- 


NOTATION FOE AERODYNAMIC LIFT AND MOMENTS 


In presenting his numerical results Dietze has introduced repre- 
sentations of lift force and moments • involving only real quantities . 
Complex notation is used in the present report in order to conform more 
nearly to current American practice; however the essential features of 
Dietze* s notation (translation of reference 3) are retained. Lift and 
moments on an airfoil strip of width Az for an airfoil with simple 
hinged flap are as follows: 
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AP S = JtpV 2 Az^CD r 2 c ss - kgs) Os + (“V^sd “ k SD) ( ?D + (“r^SR “ ^r)^ 
AMd = npV 2 (|) Az|^£D r 2 Cjjg - + (^ r 2 c DD - ^^qj) + (“r 2 c DR - k DE ^q R J- 

AM r =• *PV 2 (1) Az^aij. 2 ^ - kjg^qg + (^j. 2 c m - hRo)^ + ^ r 2 c m - *1®)^ 
where 


APg 

AMq 

AMr 



lift force, positive upward 

st alling moment on airfoil plus flap, referred to quarter- 
chord point 

hinge moment on flap, positive in same sense as AMr 
downward displacement of quarter-chord point 


SR 


rotation of airfoil in stalling direction 
rotation of flap in stalling direction 


k gh - + ^gb” 

% = (q h ‘ + iq^e 


io>t 


} g, h = S, D, R 


The quantities Cg^ may he expressed in terms of the functions 04 , 
and $22 defined in reference 6 as follows: 


C gh 



S 

D 

R 

s 

1 

1/2 

$ 4 / 2 ^ 

I) 

1/2 

3/8 

Qrf/hlt 

R 

04 / 2 * 


$ 12 /^-^ 
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DISCUSSION OF NUMERICAL RESULTS 


rv 


The coefficients kgg, kj)g, kgp, and k^g, which do not' depend on 
the ratio of flap chord to total chord t r , are presented in table III 
and in figures 1 to 8 for M = 0.7 and a range of £% from 0 to 0.7. 

The hinge-moment coefficients kpg and kgp associated with air- 
foil flapping and rotational motions are presented in table IV and 
figures 9 to 12 for M = 0.7 , reduced frequencies from to*. = 0 to 0-7, 
and ratios tr = 0.15, 0.24, 0.33, and 0.42. Coefficients P n in the 
series representation for y 

y Z -2v(p 0 cot sin 

\ 2 n=l 

are presented in table V for both flapping and rotational motions. It 
is noted that the coefficients kgg and kjg) may be computed for any 

ratio of flap chord to total chord without further iterations by 
inserting the coefficients P n for the appropriate type of airfoil 
motion into the formula 



c Eh 


" ^BhPr 2 + « + XZ P n' 


h « S, D 


The quantities Q n are expressed as functions of t r as follows: 

cos 0 = 2tjj - 1, 0 ^ 0 — it 

Qq = (« - 0)(-l + 2 cos 0) + 2 sin 0 - i sin 20 

3 1 

= 2(rt - 0) cos 0 + — s3jl 0 + £ sin 30 
q 2 a -(it - 0) - £ sin 20 + sin 40 


2 sin n0 sin (n + 2)0 


o - sln ( n ~ g ) e 

(n - l) (n - 2) " (n - l)(n + l) " (n + l)(n + 2)' 


n > 2 
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The aerodynamic coefficients associated with rotational motion of 
the flap kg E , kjj E , and kjjg are presented in table VI and figures 13 

to 18 for M = 0 . 7 * reduced frequencies ranging from 'a^ = 0 to 0 . 7 * and 
ratios of flap chord to total chord t e = 0 . 15 , 0.24* 0 . 33 , and 0.42. 

From five to eight iterations have been employed in the calculation 
of 7 . In computing the coefficient kpR the accuracy depends on the 

number of coefficients used as a starting basis for 7 ^^ . The number of 

coefficients is reduced by 3 in each successive iteration; in the calcu- 
lations described in this report 22 coefficients have been used as a 
starting basis. All coefficients available at a given stage of the 
iteration have been used to calculate the contribution of y n to y 

and the contribution of 7 inc has been obtained in closed form from the 
results of reference . 6 . 

In recalculating the coefficients which are independent of t e and 
of the remaining coefficients for t e = 0 . 15 , values have been obtained 
which differ in general by less than 1 percent from those given by 
Dietze. There are a few isolated exceptions, however. An error of 
7 percent has been found in the imaginary part of kpg for ojj. = 0 . 10 . 

Also there are errors in Dietze ’ s values of the imaginary part of k^g 

at ajj. = 0.02 for all Mach numbers tabulated, including M = 0. t i T hf» 
entry for M = 0 should read lO^kgg'* = 1.18 instead of 1.03. 

Apparently the same error has been carried through for all Mach numbers . 

A similar error occurs in the imaginary part of kjg) at (Dj, = 0.60; the 

entry for M = 0 should be lOTcgu * » = 131.8 instead of 132.8, and 

this error has been carried through for other values of Mach number as 
well. 


Chance Vought Aircraft 

Division of United Aircraft Corporation 
Dallas, Tex., July 19, 1949 
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TABLE I.- VALUES OF AED p^ FOB M = 0.? 
jotjL -* ccj* + ioi”; Pijt = p^* + iPit*’] 
































BACA TN 2213 


21 


TABLE I.- VALUES OF a± AHD 
FOB M = 0*7 - Continued 



0Q2 

010 


002* 

002” 

010 ’ 

0io’* 

0.02 

-0.00001 

-0.00001 

-0.004i6 


0.001148 

.o4 

-.00004 

-.00009 

-.oi4oi 


.00591 

.06 

-.00018 

-.00024 

-.02809 


.01331 

.08 

-.006142 

-.ooo48 

-.04561 


.02366 

.10 

-.00079 

-.00081 

-.06601 


.03696 

.20 

-.00519 

-.00360 

-.19879 


.14786 

.30 

-.01438 

-.60769 

-.36140 


.33268 

.4o 

-.02896 

-.01257 

-. 53^16 


•59143 

.50 

-.04790 

-.01796 

-.70335 


.92411 

.60 

-.07225 

-.02363 

-.85837 


1.33071 

.70 

-.10161 

-.02952 

-.99061 


1.81125 


P11 

012 


011 ' 

011* * 

012 ' 

012” 

0.02 

-0.00001 

-0. 0l<802 

0.00047 


0 

.04 

-.00005 

-.09611 

.00188 


0 

.06 

-.00016 

-.14430 

.00424 


fmm 

.08 

-.00038 

-.19262 

.00753 


SB- 

.10 

-.00075 

-.24108 

.01177 


SB 

.20 

-.00599 

-.48583 

.04706 


SE; 

.30 

-.02020 

-.73434 

.10590 


.0 

.4o 

-.04790 

-.98513 

.18826 


0 

.50 

-.09356 

-1.23571 

.29416 


,0 

.60 

-.16166 

-1.148287 

.42358 


0 

.70 


-1.72280 

.57654 


0 


o 
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TABLE I • — VALUES OF oc^ AUD 
FOR M = 0,7 - Concluded 


“r 

013 

020 


p 13* 

p 13‘* 

p 20 f 

P 20' * 

0.02 

0 

-o.ooooo4 

0.000003 


.04 

0 

-.00003 

.00002 

.00005 

.06 

0 

-.00010 

.00008 

.00014 

.08 

0 

-.00024 

.00019 

.00030 

.10 

0 

-.0001*8 

.00037 

.00053 

.20 

0 

-.00381 

.00299 

.00288 

• 30 

0 

. -.01286 

.01010 

.00713 

.40 

0 

-.03049 

.02395 

.01251 

•50 

0 

-.05956 

.04678 

.01778 

.60 

0 

-.10292 

.08083 

.02134 

.70 

0 

-.16343 

.12836 



“r 

021 

p 22 


1 

P 23 


021 ' 

021” 

022* 

022” 

023 J 

023” 

0.02 

0 

-0.01*801 

-0.00094 

0 

0 

-0.000002 

.04 

0 

-.09602 

-.00376 

0 

0 

-.00002 

.06 

0 

-.14404 

-.00847 

0 

0 

-.00005 

.08 

0 

-.19208 

-.01506 

0 

0 

-.00012 

.10 

. 0 

-.24015 

-02353 

0 

0 

-.00024 

.20 

0 

-.48101 

-.09413 

0 

0 

-.00191 

• 30 

0 

-72331 

-.21179 

0 

0 

-.00643 

.4o 

0 

--96774 

-37652 

0 

0 

-.01525 

.50 

0 

-1.21504 

-.58831 

0 

0 

-.02978 

.60 

0 

-1.46591 

-.84716 

0 

0 

-.05146 

.70 

0 

-1.72106 

-1.15308 

0 

0 

-.08171 



TABU5 H.- VALUES OF KERNEL AMD KERNEL DHTEHENCE 


B 

K«(s,0.7) 

K' '(a,0.7) 

K«(s,0) 

K’’(s,0) 

££’(8,0.7) 

AK”(b,0.7) 

£^'(8,0.7) 

£^"(8,6.7) 

£Kg’{8,0.7) 

££2' '(b»o-7) 

-2.040000 

0.011560 

0.012424 

0.015418 

-0.022375 

-O.OO3858 

0.034799 

-0.077960 

0.170m 

0.074101 

-0.135752 

-1.748971 

.001205 

.021261 

1 021173 

—027647 

-.019968 

.048909 

-099707 

,l4o884 

•079739 

-.091975 

-1.457143 

-.016623 

.022300 

.030307 

-.035040 

-.046930 

.057339 

—120836 

.109422 

.073896 

-.052083 

-1.238571 

-.032406 

.013815 

.041122 

-.042762 

-073528 

.056578 

-136632 

.084184 

.063105 

-.027606 

-1.020000 

—047369 

-.007936 

.058209 

-.053466 

-105577 

.045530 

-153120 

.056930 

.047542 

-,on4oo 

-.001429 

-051371 

-.040850 

.087487 

-069077 

-.138858 

.028227 

-.171854 

.026681 

.03^96 

.001546 

-582857 

-.034254 

-.096274 

.144397 

-093655 

-.178651 

-002619 

-.197094 

-.008492 

.018443 

.005873 

-.437143 

.006705 

-.150654 

.220011 

-.119616 

-.213307 

-.031038 

—223823 

-.03674a 

,010516 

.005704 

-.291429 

..112041 

—230339 

.381832 

-.161756 

-269791 

-.068583 

-.274424 

-.072497 

.004632 

.003913 

-.218571 

.228886 

-.288192 

.550946 

-.195116 

-.322060 

-.093076 

-.324627 

-.095786 

.002566 

.002710 

-.145714 

.474420 

-371253 

.899102 

-.246179 

-.424681 

-.125074 

-425794 

-.126578 

.001113 

.001504 

-.072857 

1.237186 

-.516904 

1.970378 

-.341703 

-.733192 

-175201 

-733459 

-.175698 

.000267 

.000497 

-.043714 

2.268691 

-.626625 

3.415745 

-.416610 

-1.147054 

-.210016 

-1.147150 

—210226 

.000096 

.000210 

.043714 

-2.968175 

-.610725 

-3.915267 

-.394760 

,947092 

-.215966' 

•94703J0 

-.216180 

.000062 

.000214 

.072857 

-i. 935504 

-.490413 

-2.469051 

-305307 

.533547 

-185106 

.533319 

-.385622 

.000228 

.000516 

.145714 

-1.167286 

-318355 

-1.393803 

-.173580 

.226517 

-144775 

.225654 

-.146425 

.000862 

.001650 

.218571 

-.912709 

-.209057 

-1.0390*50 

-.086699 

.126341 

-.122358 

.124487 

-125556 

.001854 

.003197" 

.291429 

-.783300 

-.125224 

-.860749 

-.018096 

.077449 

-.107128 

.074284 

-.112190 

.003165 

.005062 

.437143 

-.642669 

.005298 

-672994 

.092060 

.030324 

-.086762 

.023683 

— 096282 

.006642 

.009520 

.582857 

-553769 

.108392 

-.561844 

.181551 

.008075 

-.073158 

-.003046 

-.087878 

.011121 

.014720 

.801429 

-.444153 

.231165 

-435327 

.290098 

-.008826 

-.058933 

-.028286 

-.082475 

.019466 

.023542 

1.020000 

-.336260 

.321970 

-.319892 

.372588 

-.016369 

-.050635 

-.047020 

-.081996 

.030652 

.031379 

1.238571 

-.226728 

.389086 

-.204195 

.429897 

-.022532 

-.040811 

-.063508 

-.084512 

.040973 

.043702 

1.457143 

-.112563 

.427351, 

-.087011 

.461734 

-.025552 

-.034383 

-.079314 

-.089044 

.033762 

.054660 

1.748571 

.039420 

.437118 

.067247 

.464472 

—027827 

-027354 

-.100433 

-.097275 

.072605 

.069921 

2.040000 

.181763 

‘ .402045 

.210670 

.423589 

-.028906 

-021544 

-.122180 

-.107301 

.093274 

.085757 
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TABLE HI.- AERODYNAMIC COEFFICIENTS, I ND EP E ND ENT OF BATIO OF FLAP 

CHORD TO TOTAL CHORD 

« 

[m - o.fj 


Coefficient^ 

0 

0.02 

0.04 

m 

. o', op 

B 

D 


JBjt 

B 

m 

0.7 

kgg* X 10 4 

0 

75.0 

216.2 

373.8 

534.6 

688.6 

1347 

1893 

2432 

3034 

3798 

4767 

E aa M X 10 4 

klU 

0 

517.0 

952.4 

1323 

1650 

1944 

3188 

4334 

« 

5530 

6838 

8252 

978 9 

kjjg* X 10 4 


-3.0 


-23.5 

-38.6 

-56.1 

-173.3 

-337.4 

-546.2 

-791.2 

-1010 

-1173 

IqB*' X 10 4 

H 

0.5 

H 

.8,2 

16.3 

26.9 

115.3 

271.8 

522.2 

904.7 

1459 

2222 

^sd' x io 4 

VO 

o 

o 

.. 

& 

25,937 

24,025 

22,430 

21,191 

20,186 

17,461 

16,676 

16,802 

17,487 

18,569 

19,909 

*sd" X10 4 

0 

-3046 

-4058 

-4323 

-te53 

-3969 

-1659 

750.5 

2929 

4851 

6468 

7757 

* DD ' X 10 4 


25.5 

72.8 

125.6 

179.1 

236.9 

508.9 

838.0 

1301 

1972 

2883 

4136 

W X 10 4 


399-4 

770.9 

1122 

l46o 

1792 

3to2 

5044 

6760 

8532 

10,302 

n,973 
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IV.- HOTOE -MOMENT COEFTIdEWB TOR FLAPPING AND ROTATION OF AIRFOIL 
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TABLE V.- VALUES OF COEFFICIENTS P n IB SE RIES REPRESENTATION 
7- = -2V |p 0 cot ^ + 2 Pr> sin n$^ 


^0 co ^ g + 2 ^ P n sin n$^ 

[p a .p n * + iv] 



For flapping of airfoil 


aij. = 0.02 ajp = 0.04 


CDp = 0.06 



0.001*05 - 0.02580 0.01195 0.0^726 

-.00050 .00005 -.00194 .00036 


0.02102 0.06523 

-.000413 .00091 


0 

0 

0 


tttj. = 0.08 


0 

0 

0 


air = 0.10 


.00003 


.00002 


.00009 


aij. = 0.20 


a* = 


Pn’ 

Pn 1 ' ‘ 

1 

Pn" 

1 

Pn" 

Pn’ 

Pn" 

0.03053 

O.O8067 

0.03991 

0.09420 

0.08342 

0.14570 

0.12386 

0.18312 

-.00700 

.00181 

-.0101*8 

.00302 

-.03610 

.01371 

-.07419 

.03356 

.00006 

.00018 

.00013 

.00034 

.00123 

.00218 

.00454 

.00039 

0 

0 

.00011 

-.00004 

.00013 

-.00007 

.00057 

-.00043 

0 

0 

0 

0 

0 

0 

-.00004 

-.00003 


a^, = 0.40 


aip = 0.50 


(Dp = 0.60 



0.16431 

0.21075 

0.20492 

0.22786 

0.24175 

0.23334 

0.27133 

0.22743 

-.12273 

.06574 

-.17824 

.11402 

-.23186 

.17927 

-.27797 

.26203 

.01189 

.01353 

.02588 

.02355 

.04912 

.03341 

.08434 

.03984 

.00167 

-.00157 

.00378 

-.00432 

.00634 

-.01014 

.0081*8 

-.02027 

-.00018 

-.00014 

-.00059 

-.00044 

-.00165 

-.00096 

-.00355 

-.00157 
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TABLE V.- VALUES OF COEFFICIENTS P n IK SERIES REPRESENTATION 
7 = -2V (pq cot ^ + 2 P n sin - Concluded 


For rotation of airfoil 



CDj. = 0.02 


CDy = 0.04 


OJy = 0.06 



291*01 -0.19220 1.19290 -0.27984 1.10671 -0.32806 


00272 


.03991 • .00793 


.07696 .01391 


.11189 


00032 -.00003 .00125 -.00013 .00270 -.00035 



= O.30 




0 1.03781 
.02012 
.00461 

0 
0 


. 35791 * 
.11*52 9 

-.00069 

-.00009 

0 


■‘•n 


.97989 

.02691 

.OO697 

-.00001* 

0 


-0.37637 

.17790 

-.00126 

-.00018 

0 


‘n 


0.80131 

.06175 

.02586 

-.00037 

-.00001* 


-0.1*1701 

.331*05 

-.00612 

-.00135 

0 


OXp = 

= 0.50 

p 1 

P T ’ 
r n 

0.5561*1* 

-0.5l*19l* 

.2551*3 

-781*1*9 

.15191* 

-.06873 

-.01190 

-.02105 

-.00196 

-.00169 


p * 
r n 

p ’ * 
r n 

0.701*75 

-0.1*5044 

.10655 

.1*8796 

.0561*9 

-.01645 

-.00159 

-.00454 

-.00029 

.00014 

CBj. = 

= 0.70 

P * 

p * » 
r n 

0.371*84 

-0.61770 

.49812 

1.00553 

.26823 

-.19176 

-.04688 

-.05153 

-.00929 

.00848 


= 0 . 1*0 


^n 


0.63101 

.16912 

.09907 

-.001*72 

-.00097 


-0.1*9399 

.61*01*3 

-.03557 

-.01080 

.00055 


^n 


0.1*7233 

.36613 

.21279 

-.02383 

.00531 


= 0.60 


p n 


-0.59052 

.91393 

-.11630 

-.03601 

.00378 























































iakz vi.- jmwnMKM ocmnamns mm w rus fitnuiov 


i 

i 


U* - o.7j 


OoafYldsiit 

\«r 

t r\ 

0 

OrOS 

0 .0* 

0.09 

0.06 

0.1 

0.8 

0.3 

o.k 

0.9 

0.6 

0.7 

igB' x id' 

0.19 

13 , *58 

10 , *11 

11,398 

10,38* 

9901 

919* 

7380 

6373 

9180 

■3319 

*901 

**79 

kga" X 10* 


0 

-1539 

-8990 

-3183 

-3330 

-3*39 

-3*02 

-31fl7 

- yA & 

-2961 

-8917 

-SB69 

W X ID* 

O.B* 

16,7*8 

15,381 

1 *,0O6 

13,1*9 

18,867 

11,338 

9313 

8296 

7716 

731* 

6963 

66 aB 

W * 10* 


0 

“2233 

-3S09 

-3®T 

-3930 

-*083 

-3769 

-3336 

-3031 

-8889 

-8906 

-8776 ’ 

kgg' X 10* 

0.33 

19,29* 

17,919 

19,*03 

13,809 

1 *,8UB 

13, *09 

10,916 

9986 

9**J 

916* 

8960 

888 * 

Kgj" X 10* 


0 

-8*97 

-33*8 

-*0*3 

-*838 

-*316 

-36U 

-3137 

-865* 

-8337 

-8163 

-806* 

*bb' x 

0.*8 

01,370 

19,7*3 

lfl,803 

19,919 

13,9*6 

1*,9T9 

18,*** 

11, *06 

11,018 

10,909 

10,9*9 

11,020 

km" X 10* 


0 

-807* 

-3133 

-*889 

-**17 

-**83 

-368* 

-2723 

-196* 

-1**6 

-1116 

-963.6 

kja' x 10* 

0.13 

3*11 

3*83 

3*53 

3*BT 

5388 

3399 

9138 

3918 

607* 

6602 

fiseo 

6877 

<W' x 


0 

93.* 

m .9 

1*9.3 

176.0 

19*.? 

883.8 

177.3 

39.6 

-122.2 

-36*. 6 

-653.1 

km 1 x io* 

0.2* 

3797 

3908 

3639 

3978 

3988 

3997 

680* 

6*66 

6768 

7069 

73*8 

7386 

km" x io* 


0 

122 . B 

223.3 

313.1 

393.8 

*63.* 

739.3 

9*9-* 

10*9 

1063 

990.3 

609.9 

W xio* 

0.33 

3617 

393* 

3*73 

3721 

3767 

3983 

6100 

6*19 

6838 

788 * 

7769 

8262 

km" x io* 


0 

lftL.9 

3*8.* 

*97-3 

630.1 

7**.0 

186* 

176* 

2160 

8*93 

8633 

8686 

Wxio* 

O.kfi 

310 * 

3U3 

3169 

3806 

3871 

3338 

36*0 

6009 

6*66 

7038 

7708 

8**7 

kQ," x io* 


0 

837.9 

*38.9 

650.3 

936.3 

m3 

1819 

2362 

3873 

3878 

*373 

*782 

km 1 x io* 

0.13 

161.7 

179.T 

179.1 

17*. 0 

178.* 

171.1 

l£8.* 

163.6 

170.7 

173.8 

176.O 

178.9 

km" x ID* 


0 

-2.X 

-out 

Q,k 

8.7 

3.1 

19.0 

32.1 

**.3 

33.9 

66.2 

73.9 

k^-xio* 

0 .2* 

Me . 3 

*70.8 

* 6 oa 

*31.3 

***.6 

*39.3 

*87.7 

*29-3 

*36.9 

*98.0 

*67.3 

*e*,2 

*m" Xio* 


0 

-9.0 

-9,8 

0.9 

10.0 

20.1 

73-0 

130.* 

181.6 

228.6 

272.3 

318.7 

kjg,' x ID* 

0.33 

930.* 

91 9.9 


868 .* 

830.3 

936,7 

903.3 

809.9 

836.3 

978.9 

919.6 

9T*.5 

km" x io* 


0 

-8*.8 


0.2 

23.2 

*9.0 

,169.8 

389.6 

*66.1 

991.9 

708.3 

917.1 

km' xio* 


l£08 

13*3 

1*83 

1*36 

1*01 

V73 

1311 

1319 

1361 

1**8 

1393 

168* 

km" xio* 


0 

-31.7 

- 59 a 

-4*1 

tea 

9*. 3 

376a 

668.8 

9*8.2 ■ 

1206 

1 *** 

1690 


ro 

03 


t 
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Figure 1.- Real part of kg g against reduced frequency for . M = 0.7. 
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Figure 5.- Real part of for M = 0.7. 
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Figure 7.- Real part of k DD for M = 0.7 


36 


MCA TN 2213 



Figure 8.- Imaginary part of 


for M = 0.7, 
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Figure 10.- Imaginary part of for M = 0.7 



Figure 11.- Real part of for M 
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for M = 0.7, 
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Figure 16.- Imaginary part of kp R for M = 0.7. 







